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We present a model which predicts inflation without the presence of inflaton fields, based on the 
eR^ and Starobinsky models. It links the above models to the reheating epoch with conformally 
conpled massive particles created at the end of inflation. In the original Starobinsky model, the 
reheating era was created by massless non-conformally coupled particles. We assume here that 
non-conformal coupling to gravitation does not exist. In the tP? model, inflation is produced by 
the gravitational Lagrangian to which a term eR^ is added, where e is a constant and R is the Ricci 
scalar. Inflation is created by vacuum fluctuations in the Starobinsky model. Both models have 
the same late-inflation time-dependence, which is described by a characteristic mass M. There is 
a free parameter Ho on the order of the Planck mass Mpi that determines the Hubble parameter 
near the Planck epoch and which depends upon the number and type of particles creating the 
vacuum fluctuations in the Starobinsky model. In our model, we assume the existence of particles 
with a mass m, on the order of M, conformally coupled to gravity, that have a long decay time. 
Taking m = FM, we investigate values of A = 0.5 and 0.3. These particles, produced ^ 60 e-folds 
before the end of inflation, created the nearly scale invariant scalar density fluctuations which are 
observed. Gravitational waves (tensor fluctuations) were also produced at this epoch. At tend, the 
Hubble parameter begins to oscillate rapidly, gravitationally producing the bulk of the m particles, 
which we identify as the origin of the matter content of the Universe today. The time required for 
the Universe to dissipate its vacuum energy into m particles is found to be tdm — &Mli/M^F. We 
assume that the reheating time Irh needed for the m particles to decay into relativistic particles, is 
very much greater than that necessary to create the m particles, tdis- A particle physics theory of m 
can, in principle, predict their decay rate Fmr = tpn- From the ratio / = tdis/tuH, F and g, (the 
total number of degrees of freedom of the relativistic particles) we can, then, evaluate the maximum 
temperature of the Universe Tmax and the reheat temperature Thh at Irh- From the observed 
scalar fluctuations at large scales, 5p/p ^ 10~®, we have the prediction M = 1.15 x 10~®Mpi and 
the ratio of the tensor to scalar fluctuations, r = 6.8 x 10“"^. Thus our model predicts M, tdis, tend, 
Tmax, Trh, tmax, and Irh as a function of /, F, and p* (and to a weaker extent the particle content 
of the vacuum near the Planck epoch). A measured value of r that is appreciably different from 
r = 6.8 X 10“"^ would discard our model (as well as the Starobinsky and eR^ models). 

PACS numbers: 98.80.-k, 98.80.Es 


I. INTRODUCTION 

The standard model of inflation, based on the exis¬ 
tence of a scalar inflaton held, makes the following as¬ 
sumptions: 
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1. The beginning of inflation occurs at an energy << 
Mpi (Planck energy). Its origin is unknown and 
the state of the Universe before the beginning of 
inflation is undehned; 

2. A large initial displacement of (j) from the minimum 
of ¥{(/)) is necessary for the onset of inflation (such 
as in the chaotic inflation model); 

3. The potential energy of the inhaton dominates its 
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kinetic energy during inflation; and 

4. The inflaton potential, V((/>), and its first derivative 
are defined by observations ~ 60 e-folds before the 
end of inflation. 

In complex inflation theories, there can be more than one 
inflaton. (See P| for a recent review of inflation theory 
with inflatons.) 

Here we present a model which links the Starobinsky 
and models of inflation, where R is the Ricci scalar 
to the reheating era. All three models, the Starobinsky, 
and ours, do not involve inflatons to create inflation. 
They also avoid most of the above assumptions. 

In the Starobinsky model Q, an term in the ef¬ 
fective L^rangian dominates inflation at late times (see 
also 1^, P]). There is no sharp boundary between the 
Starobinsky model and model since the latter is the 
particular case of the former in the limit M <C i?o (us¬ 
ing the notation of Eq.(9), with some small non-local 
terms (due to non-zero rest masses of conformally cou¬ 
pled quantum fields) omitted. However, the Starobin¬ 
sky and eB? models have the same qualitative behav¬ 
ior at ~ 60 e-folds before the end of inflation, when the 
presently observed scalar and tensor fluctuations were 
produced. 

Both the Starobinsky and eB? cosmologies are char¬ 
acterized by a single mass M (= Mpi/-^2A e in the eB^ 
cosmology and Mpi/^/A&Trki in the Starobinsky model, 
where Mpi = is the Planck mass and fci is the 

coefficient of the term which contains the second deriva¬ 
tive of B in the quantum corrected vacuum expectation 
value of the energy-momentum tensor [Eq.Q]). The 
mass M characterizes the end of the inflation period, 
during which, the Hubble parameter varies slowly. A 
period then begins, in which H oscillates rapidly as 
H oc (1/t) cos^ ojt and the cosmological scale factor varies 
as a{t) oc t^/^[l-|-sin [2u;t]/(Swt)], where w ~ M/2. When 
averaged over several oscillations, the Universe expands 
as a classical matter-dominated Universe. 

Although the eB? model can be considered to be the 
simplest way to produce inflation, i.e., by means of a sim¬ 
ple modification of the gravitational Lagrangian, we con¬ 
centrate here on the Starobinsky model since it is more 


complete. It links the beginning of the inflation period 
to the beginning of the Universe and also describes the 
end of inflation in detail. 


The Starobinsky model suggests that for energy densi¬ 
ties and curvatures near the Planck scale, quantum cor¬ 
rections to Einstein’s equations become important (as 


discussed in detail by Vilenkin |3). In the Starobinsky 
model, inflation is driven W one-loop corrections due 


to quantized matter fields [3 (see also |6L iTl I8L I9L Ilf 
If~i1|b The model is consistent with a Universe that was 
spontaneously created, as discussed by Grishchuk and 
Zel’dovich |1^ . 

The beginning of the Starobinsky inflation period can 
be associated with the beginning of the Universe due to 
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quantum fluctuations of the vacuum. Tryon 
the first to suggest that a closed Universe can be cre¬ 
ated spontaneously as a result of a quantum fluctua¬ 


tion. Vilenkin 

H, 


14 ll,’i| . Zel’dovich and Starobinsky 
and Linde |17| were the first to attempt to describe 
the quantum creation of a Universe in the framework of 
quantum gravity. The picture that emerges is one of a 
Universe tunneling quantum mechanically to a de Sitter 
space time. At the moment of nucleation (t = 0), the 
Universe has a size a(0) = . This is the beginning of 

time and, from that point on, the Universe evolves along 
the lines of the inflation scenario. 


In the Starobinsky model, inflation is produced by the 
vacuum energy pv, which has negative pressure, P = 
—pv- Inflation in both the Starobinsky and our models 
can be described by an effective geometric scalar particle 
M. In our model, there is an additional massive particle 
m produced at the end of inflation, which is freely moving 
and which produces positive pressure. 


Structure in the Universe primarily comes from al¬ 
most scale-invariant superhorizon curvature perturba¬ 
tions Q, . In our model, a mass m is much less than 
the Hubble parameter during inflation. The mechanism 
of m particle production from inflation is based on the 
observation that particles that are massive in the present- 
day vacuum, could have been very fight during inflation. 
This implies that fluctuations of a generic scalar field y 
with mass m <C iL during inflation are copio usl'v 
ated, with an almost scale invariant spectrum |2C 
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The particles become heavy and non-relativistic at the 
end of inflation. 

The end of the Starobinsky inflation period has been 
suggested to be due to the masses of the particles in the 
vacuum fluctuations [14], Thus, since the mass M de¬ 
scribes the end of inflation in both the Starobinsky and 
our models, M is a natural mass scale for the particles 
that are created at the end of inflation. In our model, we 
then have the scenario that particles of mass compara¬ 
ble to the mass M in the vacuum fluctuations first create 
the inflation, after which, particles of mass m comparable 
to, but slightly less than M, are produced from the vac¬ 
uum due to the rapid change of the Hubble parameter. 
The particles m are conformally coupled to gravitation 
(Ricci scalar). These massive particles create the reheat¬ 
ing epoch of the Universe. 

Our model can be compared with that of the Starobin¬ 
sky model, in which massless particles, non-conformally 
coupled to gravitation, directly create the reheating era. 
Here we assume that non-conformal coupling does not 
exist. 

Previously, gravitational production of massive parti¬ 
cles has been investigated in order to explain the ob¬ 
served ultra-high energy cosmic rays, produced as a re¬ 
sult of heavy particle decay (masses > 10^^ GeV) IS 
The gravitational particle production of the heavy par¬ 
ticles m can be described assuming a given background 
metric jSIS- 

The paper is organized as follows. In section H, we give 
the main results of the Starobinsky model, as discussed 
by Vilenkin P] . We discuss the gravitational production 
of the m particles in section HI. In section IV, we derive 
the scalar density fluctuations produced ^ 60 e-folds be¬ 
fore the end of inflation. We obtain the ratio r of the 
tensor to scalar fluctuations in section V. The reheating 
of the Universe is discussed in section VI. Our conclusions 
and discussion are presented in section VH. 


II. THE STAROBINSKY INFLATIONARY 
MODEL 

In this section, we discuss the Starobinsky model 
of inflation, following the description and notation of 


Vilenkin. This model is based on the semiclassical Ein¬ 
stein equations, 

R^J.v - = -SttG < Tfj,^ > , (1) 

which assume a spatially flat Robertson-Walker metric, 

= a^(?7)(d?7^ — dx^), (2) 

where dt = adrj, t{rj) is the proper (conformal) time and 
a is the scale factor. 

The quantum corrections to the expectation value of 
the energy-momentum tensor for massless particles in 
curved space time are 

+ (3) 

where 

= 2R.^.„- 2gf,^ R-a'"" + , 


— Rfj,'^Rua — -^RRfiu — -^g^ivR'^^Rar + R^ ^ 

( 4 ) 

and ki , k^ are constants. The value of the coefficient ki 
can be determined by observations. However, the value 
of fca is fixed by the condition 


where iVp, iVj^/ 2 , and Ni are the numbers of quantum 
matter fields with spins 0, 1/2, and 1, respectively. 

It is convenient to define two new parameters, Hq and 
M, in terms of ki, fcp, and the Planck mass, Mpi: 


and 


Ho = 


Mpi 

VSnko 


M = 


Mpi 

>/487r/ci 


( 6 ) 

( 7 ) 


The Planck constant h and the speed of light c are given 
in natural units, h = c = 1, and Mpi = = 1.22 x 

10^®GeV. To evaluate idp, a minimal SU(5) model, where 
No = 34, Ni /2 = 45, Ni = 24, is assumed, giving 


Ho = 0.74 Mpi. 


( 8 ) 
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From the above, the energy-momentum tensor in 
Eq.Q is 


>=i 


(MelY+Wh (MelY 


(9) 

Whereas = Q for classical conformally invariant fields, 
a trace anomaly arises from the process of regularization 
(see HElEi), 

2 


< rr >= 


SttHI 


) R,,-’ 


( 10 ) 

The evolution equation for the Hubble parameter H{t){= 
a/a) in a flat Universe is 

-Hi) = {2H ii+ 6H‘^ H - H"^). (11) 

Assuming that H is slowly varying during inflation, H <C 
H^, and that H H H, the solution of Eq. cu is 


H = Hn tanh 


7- 


GHn 


where 


7 = 1/2 ln[2/,5o], 


So = \Hin - Ho\/Ho , 


( 12 ) 


(13) 


(14) 


and i7in is the initial value of the Hubble parameter in 
the inflation era. 

We see from Eq. (cai that H changes on a time scale. 


r ~ 6 Ho/M^ . 


(15) 


A long period of inflation occurs when M'^ ^ 6 Hi. The 
solution of Eq. in terms of a, valid up to the time at 
the end of inflation, when H ~ M, is 

6 7 Ho 


^end — 


M2 


IS 


a{t) = Hp 


_i/ cosh 7 \Hot 


(16) 


(17) 


vcosh [7 — t/r] 7 

where r is given by Eq. 113). 

Assuming that H Ho near the end of inflation, 
En. dl 1 1I simplifies to 

2HH + 6H^H - H^ + M^H^ = 0 , (18) 


which has the solution 


7?(i) = ^ cos^ 


Mt 


L 2 J 




and 


a{t) oc 


1-f 


3Mti 


sin [Mt]. 


( 20 ) 


Thus, from Eqs. (El) and (EUl) . H{t) and a{t) are in an 
oscillating phase at the end of inflation. At a time to ^ 
M~^, the period of oscillation is much shorter than the 
average Hubble time, 2/3 to- For a time interval to ^ 
At ^ M~^, we can neglect the power law expansion in 
Eq.(E|), so that 




( 21 ) 


In the Starobinsky model, massless particles are grav¬ 
itationally produced at the end of inflation. For a scalar 
field of mass m that satisfies the equation 


D(j) + (m^ + iR)4> = 0 , 


( 22 ) 


the field is conformally coupled if ^ = 1/6 and non- 
conformally coupled if ^ ^ 1/6. For m > M/2, the par¬ 
ticle production is exponentially depressed. In a confor¬ 
mally flat spacetime, massless conformally coupled par¬ 
ticles cannot be gravitationally produced [6|, |29|. There¬ 
fore, these particles must be non-conformally coupled. 

The oscillation term in Eq. ED is small and can be con¬ 
sidered to be a perturbation in the calculation of parti¬ 
cle production. A perturbative technique for calculating 
the production of very low mass m (<j; M ) pa rticles was 
developed by Z ePdovich and Starobinsky [30j and Birrel 
and Davies [33, treating | ^ — 1/6| as a very small param¬ 
eter. They showed that for a massless non-conformally 
coupled field (m = 0, ^ ^ 1/6), the particle production 
rate is 




1\2 

6. 


where 


/4M\ 

i? « 6 a = — (-) sin [Mt]. 

\ to / 


(23) 


(24) 


Taking the average of the particle production over the 
period of oscillation and using the fact that the particles 
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are produced in pairs with energy to/ 2 per particle, the 
average rate of energy loss is (to = FM) 

^ FM ^ FM^/^ 1\2 


and the functions Xk(t) satisfy the equation for the field 

Q. 






(25) 


Xk +k^Xk + 




a Xk = 0 . 


(32) 


The rate at which the vacuum energy is dissipated into 
particles is 


Linearizing Eq. and using Eqs. m and (P|l . we have 


T = pIp=- 


3 FM^ / ^ 1 \ 2 


2 Ml, 


(«-i)' 


where 


P = 


Q-kG 


(26) 


(27) 


Xk + Wfc Xk- 
where 


3Mto 


TO^ sin [Mt] Xfc = 0 , 


3to ’ 

(28) 

If the main contribution to the particle production comes 
from the modes with k ^ M/2 and the mass to <C M/2, 
we can replace with in Eg. 13411. The expression 


(29) 

of the production of to particles for ^ = 1/6 (conformal 
production), 


is the energy density of the particles. 
Since a oc , we have 


so that Eq. m becomes 
P = 

(the Friedmann relation). 


III. GRAVITATIONAL PRODUCTION OF 
PARTICLES 

As discussed in section II and in detail by Vilenkin, 
massless particles (<C M) are non-conformally produced 
at the end of inflation in the Starobinsky model . These 
massless particles are assumed to reheat the Universe. 

Although the main purpose of the analysis of Vilenkin 
was to evaluate the gravitational production of mass¬ 
less non-conformally coupled particles in the Starobinsky 
model [Ea. ll^ ]. an expression for the conformally cou¬ 
pled gravitational production of massive to particles was 
also obtained. We use this expression in our analysis of 
massive to particles. 

The created (p field [Eq. (I22|l ] can be expanded in terms 
of creation and annihilation operators. 


and 


Wfc = {k^ + 


to2 = to2 - 3 ( C - ^ ) m2 . 


(33) 


(34) 


(35) 


r = 


Gto^ 

6M 


(36) 


was presented in the original paper 

As noted by Vilenkin, although Eq. (EH) was derived 
assuming to <C M, it also gives a correct order of magni¬ 
tude for the conformally coupled decay rate for to ~ M, 
{m = FM, E = 0.3, 0.5) 

m3 


r = 


6 M2/ 


(37) 


We can compare the time it takes to dissipate the vac¬ 
uum energy due to the gravitational production of mass¬ 
less particles, non-conformally coupled, using Eq. 


, p-15 _ 2 M|,, 1 

J- j — g ■ 


3 F3M3 _ 1/6) 


2 ’ 


(38) 


with that of massive particles M, conformally coupled, 
from Eq. E3, 


(j){x) = d^k akUk{x) + alul{x) 


where 


Ukix) = 


(27r)3/2 




'Xk{t) 


(30) 


(31) 


f-dis — 


6 Mpi 


M3 ■ 

The ratio of the two times is 

^dis^ _ F 


(39) 


f-dis 9 (^ — 1/6) 


(40) 
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Since it is generally assumed that (^ — 1/6) <C 1, the vac¬ 
uum energy is dissipated more rapidly in the case of the 
emission of the massive m conformally coupled particles 
[Eo. iBIHi ]. 

In general, the vacuum can lose energy, both by the 
gravitational production of m particles that are confor¬ 
mally coupled or by non-conformal gravitational produc¬ 
tion of massless relativistic particles. For simplicity, we 
assume that ^ = 1/6 and, thus, that only conformal pro¬ 
duction of m particles exist. From Fqs. and iTTUI) 
(below), we have 

^dis ^ ^ 6.8 X 10^^ tpi, (41) 

where r is the ratio of tensor to scalar fluctuations and 
tpi = Mpl = 5.39 X lO^'^^sec is the Planck time. 

In this paper, we assume the possible existence of el¬ 
ementary particles of mass m = FM, which are confor¬ 
mally produced during or at the end of inflation. From 
their lifetime, we obtain the maximum temperature of 
the Universe Tmax the reheat temperature TpH and their 
respective times, tmax, and tpp. From the gravitational 
production of the m particles at ^ 60 e-folds before the 
end of inflation, we obtain the scalar density fluctuations 
(section IV). 


IV. SCALAR DENSITY FLUCTUATIONS 

Structure in the Universe primarily comes from nearly 
scale-invariant superhorizon curvature perturbations. 
These perturbations originate from the vacuum fluctu¬ 
ations during the nearly exponential inflation. In our 
model, structure is due to a scalar field x of mass 
m <C Hqo, where Hqq is the Hubble parameter ~ 60 e- 
folds before the end of inflation. The scalar field x of 
mass m Hqq in a quasi-de Sitter phase was shown to 
produce quantum fluctuations, whose power spectrum is 
scale invariant if they are superhorizon j2ft |2l| . 

During the inflationary epoch, the fluctuations of the 
field X of mass m obey the equation 


6xk + SHSxk + ( “ ) ^Xk = 0 


(42) 




191. The fluctuations Sx are described in terms of the 


variance. 


< X >= 


d^k 

(27r)3 




(43) 


which obeys the equation 

d < x^ > 


dt 


47r^ 


(44) 


during inflation 


dependently obtained in |36l l3 


[^ . (The formula Eo. iEI) was first 


H^Hn 


^.) In our model, 
MH 


7- 


6i?o 


(45) 


until the end of inflation, before oscillations set in. Sub¬ 
stituting Eq. 103) into Eq. 031, we obtain 

3 


d < x^ > 


Hh^ 


1 - 


t 

^end 


(46) 


dt 

where tend is given by Eq. 03- 

We note that the major contribution to the variance 
comes from t tend, while very little comes from 
t ~ tend- Integrating Ea. lBHIl from (l/60)tend to tend, 
we estimate the variance of the fluctuations that were 
created ^ 60 e-folds before the end of inflation, 


2 3 

87r2 M2 

From Eq. 011, the Hubble parameter is 
-H^end — 


3 ^end 


(47) 


(48) 


at the end of inflation. Substituting Eq. llTT)ll into Eq. 0H1, 
we find 


2 M^ 

Hend — 77" 


9 xHq 


In a flat Universe, we have 
8tt 


= 


SMpi 


[PUnd + PMend] > 


(49) 


(50) 


where pv^nd the vacuum energy density at the end of 
inflation. From Eqs.011 ISOl, PVend i® given by 


PUnd = 


1 M^Mj.; 


(51) 


547r ’ 

assuming that is very much greater than pM at the 
end of inflation. 

The long wavelength x modes satisfy the equation 


5xk + 8H5xk + m^Sxk = 0 


(52) 
















and that of the inhomogeneous form, 
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Q, 


which gives 


6xk oc a cos mt. 


(53) 


They are non-relativistic, behaving like a classical homo¬ 
geneous field. Their number density is given by 


dmend rn<x > 


m 


From Eq. 63 , we then have 


Pmend ~ 




(54) 


(55) 


We assume that the m particles decay into relativistic 
particles with a decay rate, F^r Fvm = ^dis- Thus, 
we can separate the decay of the vacuum into m particles 
(Tvm) from the decay into relativistic particles. 

For the production of the m particles, we have the 
relation 


Pm T Pm — PV 7 


(56) 


where the second term on the left describes the dilution 
of the m particles due to the expansion of the Universe 
and the term on the right, the production of m particles 
due to the decay of the vacuum. 

The production of m particles starts at the end of the 
inflation period, when 


and 


py(i = tend) = Pi4nd (57) 


Pm(t — tend) — Pme„d ^ PVe„d ■ (58) 


During the time m~^ t tg, when the m par¬ 
ticles were produced, the average Flubble value was 
H = 2/3<o- The vacuum energy decays in a time 
tdis = ^vrny during which, the m particles are produced: 
pv{t) = exp [-Fy^t]. Let us take to - tdis- We 
can then describe the production of the m particles by 
Ea. (l^ in the form 

Pm T 3iFpm = P\4nd/tdis ■ (59) 

The solution of the homogeneous form of Eq. is 

— —tend) (60) 


Pm{t) = . (61) 

^dis 

We take the initial time for the decay of the m particles 
into relativistic particles to be 


^mri — ^dis ^end ■ 

From Eq. EJ, we have 

Pm{i — ^mri) = Pmi — PVend ' 


(62) 


(63) 


Since the Universe had an average cosmic scale fac¬ 
tor a oc from tend to tdis, it expanded by a factor 
(tdis/tend)^/^ and diluted pmend by a factor (tend/tdis)^ 
during this time interval. Observations show that the ra¬ 
tio of the energy density of the mass fluctuations at tdis, 
Prriend (^snd/tdis) , tO that of matter, Pmil 70 

From Mukhanov and Chibisov (1981) we have an 
approximately flat scalar fluctuation spectrum for the 
Starobinsky model with an amplitude 

1/2 


(t) 

From Eqs.(62) and (63), we then have 

(^end/^d is) Pmend 


PVend 


~ 10 ~ 5st ■ 


(64) 


(65) 


Although 5st was evaluated for massless particles, sim¬ 
ilarly to Vilenkin [8], as discussed in Sec.II, we assume 
that this is also the amplitude for the production of mas¬ 
sive particles on the order of M. From Eqs.(64) and (65), 


M ~ 1.15 X 10“^Mpi. 

V. RATIO OF TENSOR TO SCALAR 
FLUCTUATIONS 

From Vilenkin, the tensor power spectrum is 
PT{k) = \hu\^ = ^ 


( 66 ) 


27r^ 


(67) 


and the scalar power spectrum is 


Ps{k) = 5i = 


2 _ {^P/p\or 


k^ 


Pm{t) — Pmend® 


( 68 ) 










where k is the wavenumber of the fluctuations and 
{6p/ p)hor the density fluctuations on the order of the 
horizon. The ratio r of tensor to scalar fluctuations is, 
then, 

^ = 1 (69) 

- Psik) ^Sp/p)l^ ■ 

We used the observed value of the scalar fluctuations, 
(i5p/p)jjQj. ~ 10“®, to obtain the value of M in terms of 
r, 


M ~ 4.4 X Mpi. (70) 

From Eq. we obtain tdis- We find that for F = 0.5, 
tdis — 6.28 X lO^^tpi and that for F = 0.3, tdis — 4.85 x 
10^ Vi. 

We have 


r ~ 6.8 X 10"^ 


(71) 


from Eqs. fTnji and 

A small value for r, bounded from below by r > 
3 X 10“® (unless V jV in the inflaton potential is un¬ 
reasonably large), was previously indicated in [^ . 

From Eqs. (EHl and ® , PKnd and [Eqs.(|nil and 
(1^ respectively], tend [Eq.l^B))], and tdis lEa. lT?!^ ]. we 
obtain 


6.7 X lO-V-5/4 

Ho/Mpi 


(72) 


From Eqs. IZ 3 , (IISIi, and Cl for the SU(5) model, the 
initial value of the Hubble parameter Hin of the inflation 
era for F = 0.5 and 0.3 is 


V Mpi. (73) 


VI. REHEATING OF THE UNIVERSE 


In this section, the reheating of the Universe in our 
model is discussed. It is based on the discussion of this 
epoch in 0 

The particles m have a decay rate into relativistic par¬ 
ticles, 


F 


mr 


f^Vm, 


(74) 


where we assume that 


/«!. 


( 75 ) 


The equation describing the decay of the m particles 
is 


Pm T Pm — PmXmr t 
which has the analytical solution, 

Pm=Pmi(^)^ ^ ( 77 ) 

where aj is the cosmic scale factor at t = tmri- For the 
radiation energy density, pr, the evolution equation is 

Pr F Pj- = PmFmr , 


where pr is the energy density of the relativistic decay 
products. In order to obtain the energy density of the 
relativistic particles as a function of time, we also need 
the Friedmann equation. 


H^ = ^{Pm+Pr). (79) 

From t = tmri{= tdis) until t = tpp(= F”),), the m 
particles dominate the mass density and the Universe 
is matter dominated, a(t) oc with ~ 0 at tmri- 
During the m-dominated epoch, an approximate solution 
for Pr is given by 


Pr{t) 


Pmi O^mrtmri) 
lOnt 


1 - 



-5/3" 
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Pmi (Fmr^mri) I 


ai 


-3/2 


1 - ( — 
ai 


-5/2' 


(80) 

Thus, Pr rapidly increases from ~ 0 at tmri to a 
value of ~ Pmi (Fmrtmri) at tpH, decreasing thereafter 


0 


as a 


-3/2 


Once the relativistic decay products interact suffi¬ 
ciently, they thermalize and we have 


Pr = 5*7r^r'^/30, 


(81) 


where p* is the number of relativistic degrees of freedom 
and is generally estimated to be 100 F < 1000. What 
is commonly called the reheat temperature, Tpp, is not 
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FIG. 1: a) Top: The evolution of the radiation energy density pr fEo. 1181)11 1 (in units of Mpi) as a function of time (in units of 
fpi) for / = 10“®, starting from tmri (when pr — 0) for F = 0.5 (solid curve) and F = 0.3 (dashed curve), b) Botton: The 
temperature (in units of g^^^^GeV) dependence on time fEas. 18111 and ®] during the decay of the m particles into relativistic 
particles. 


TABLE I: The maximum Tmax and reheating Tpu temperatures (in units of ^y^^GeV), the times tmr-i, tmax, and tnH, for 
/ = 10~® and 10~^° and F = 0.5 and 0.3, where is the number of relativistic degrees of freedom. 


F 

/ 

Tm.. [gV^GeV] 

Trh [s^^GeV] 

tmri 

^meix [^Pz] 

tRH [tpl] 

0.5 

10-° 

~ 1.14 X 10^2 

~ 5.93 X 10“ 

~ 6.28 X 10“ 

~ 1.13 X 10“ 

~ 4.85 X 10“ 

0.5 

o 

1 

O 

~ 6.39 X 10“ 

~ 1.39 X 10® 

~ 6.28 X 10“ 

~ 1.13 X 10“ 

~ 4.85 X 10“ 

0.3 

10-° 

~ 8.26 X 1044 

~ 3.19 X 10“ 

~ 4.85 X 10^° 

~ 8.74 X 10“ 

~ 6.28 X 10“ 

0.3 

o 

1 

O 

~ 4.64 X 10“ 

~ 1.01 X 10® 

~ 4.85 X 10“ 

~ 8.74 X 10“ 

~ 6.28 X 10“ 


the maximum temperature of the Universe, T^ax, which 
is given by 

Thiax = 0.8 ( PmiTmrtmriy^'^ 9*^^^ ■ (82) 

The reheat temperature at the beginning of the radiation 
dominated epoch, = T{tpiH = r“(,), is given by 

Trh 0.55 ( (83) 

and the ratio of T^ax to Trh is 

Tmax il^Vm ^ 1.45 ^ 

since Tmr = f^Vm, where / <C 1. In Tabled we show 
the values of T^ax [Ea. ljH^ ] and Trh [Eq.(|HSI)] for r = 
6.8 X 10-4, / = 10“^ and 10-4°, and F = m/M = 0.5 
and 0.3. 

In Fig. 1, we show pr [Ea. dHIHl ] and the temperature 
as a function of time fEo. 18111 ] for / = 10-° and F = 0.5 
and 0.3 during the time interval tmri < t < tRH- The 
times tmri, tRH, and tmax are shown in Tabled 


VII. CONCLUSIONS AND DISCUSSIONS 

We presented here a model which relates the Starobin- 
sky and eR^ models, both of which predict inflation, to 
the reheating era by a massive conformally coupled par¬ 
ticle. In the original Starobinsky model, the coupling to 
the reheating era was made by massless non-conformally 
coupled particles. Here we assumed that non-conformally 
coupling to gravitation does not exist. 

In the Starobinsky model, inflation is due to vacuum 
fluctuations. Inflation is due to a modification of the 
gravitational Lagrangian in the eR? model. In both mod¬ 
els, inflation takes place without the need for a scalar 
inflaton field. 

The end of the inflation period predicted by the 
Starobinsky and tR? models is characterized by a param¬ 
eter M (M = Mpi/^/ASirki in the Starobinsky model, 
where ki is the coefficient of the tensor (EqsElll of 
the expectation value of the vacuum energy-momentum 
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tensor due to vacuum fluctuations and M = Mpij^f^t 
in the eB? model). 

In our model, we assumed the existence of particles 
of mass m and linked their conformal gravitational pro¬ 
duction during inflation to the observed scale invariant 
density fluctuations. Their gravitational production at 
the end of inflation was linked to the matter density of 
the Universe. 

Our model of two massive scalar particles M and m 
has a certain resemblance to the double inflation model 
of Gottlober, Muller, and Starobinsky [^. Gottlober, 
Miicket and Starobinsky discussed the confrontation of 
the model with observations In the simple chaotic 
inflation model, the scalar potential, (I/2)m^0^, is char¬ 
acterized by the mass m. The gravitational term, 
is characterized by the mass M in the late in¬ 
flation in the Starobinsky model. In the double inflation 
model of Gottlober, Muller and Starobinsky, both the 
scalar field potential, and the gravitational 

term, are present. Thus this double inflation 

model is described by two masses M and to. It is to be 
noted that in our model, described by the two masses 
M and to, only the mass M describes the inflation era, 
whereas the mass to is a free particle, produced at the 
end of inflation, linking the inflation era to the reheating 
era. The to particles in our model, have a decay time into 
relativistic particles tnu very much greater than tdis, the 
time for the production of the to particles at the end of 
inflation. 

Our model depends upon the parameter Hq, the nor¬ 
malization of the Starobinsky inflation solution for the 
Hubble parameter near the Planck era fEg. 11211 ]. The 
parameter Hq is on the order of Mpi and depends on 


A^Oj Ni/ 2 , and the number of quantum matter fields 
in the vacuum of spin 0, 1/2, and 1, respectively. To 
evaluate ifo, we assume a minimal SU(5) particle content 
and obtain Hq = Q.lAMpi [Eq.®]. The time at the end 
of inflation depends on Hq and M. For r ~ 6.8 x 10“^ 
and M ~ 1.15 x 10“® Mpi ~ 1.4 x 10^^ GeV, as pre¬ 
dicted by the Starobinsky model for F = 0.3, for ex¬ 
ample, we find tend — 1-37 x 10^^ tp; ~ 10“^^ sec and 
tdis — 4.8 X 10^° tpi « 10“^^ sec, respectively. 

In the future, the particle to could be incorporated into 
a particle physics theory that would define the decay time 
into relativistic particles tpp, the reheating time. From 
/ = tdis/tpH, we evaluate the maximum temperature of 
the Universe Umax and the reheat temperature Tpp as a 
function of (the number of degrees of freedom of the 
relativistic particles) and F = m/M. The times tMri, 
imax) and tpp can, then, also be evaluated. A measured 
value of r, the ratio of tensor to scalar fluctuations, that 
is appreciably different from r = 6.8 x 10“'* [Eq.(72)] 
would discard our model (as well as the Starobinsky and 
models). 

In order not to overproduce gravitinos, it is frequently 
suggested that Tpp ^ 10® GeV 1^. If this is, indeed, a 
true upper limit for Tpp, it puts limits on the possible 
values for / and F from the above analysis. However, 
it is to be noted that, for our model, 10® GeV is not a 
strong upper limit for Tpp since supersymmetry is still 
not a well developed theory. 
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